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Abstract 

The possibility of nonlinear pulses generation in Yttrium Iron Garnet thin films for arbitrary 
direction between waveguide and applied static in-plane magnetic field is considered. Up to now 
only the cases of in-plane magnetic fields either perpendicular or parallel to waveguide direction 
have been studied both experimentally and theoretically. In the present paper it is shown that also 
for other angles (besides or 90 degrees) between a waveguide and static in-plane magnetic field 
the stable bright or dark (depending on magnitude of magnetic field) solitons could be created. 

PACS numbers: 85.70.Ge; 75.30Ds; 76.50.+g 
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I. INTRODUCTION 



The investigation of magnetostatic envelope solitons in yttrium-iron garnet thin magnetic 
films is one of the "hot topics" in nowadays physics. The advanced instrumentation for mi- 
crowave pulse generation, detection and analysis together with the solid theoretical base has 
led to a growing interest in studying such localized objects. The definition "magnetostatic 
soliton" refers to a propagating pulse formed by large wavelength spin excitations which do 
not "feel" the exchange interaction and only the dipolar interactions could be taken into 
account. Thus the processes are characterized by the Landau-Lifshitz and magnetostatic 
equations. 

The linearized solutions of these equations were obtained by Damon and EshbachB 40 
years ago for arbitrary direction between wave vector of spin excitations and in-plane mag- 
netic field. The nature of those excitations has been also studied experimentally!. The 
weakly nonlinear limit for the mentioned equations also was considered for the particular 
cases when wave vector of spin excitations is either parallel (backward volume waves) or 
perpendicular (surface waves) to the direction of in-plane magnetic field. It was found out! 
that the envelope of spin excitations in both cases satisfy 2D Nonlinear Schrodinger (NLS) 
equation which permits well-known ID soliton solutions depending on the relative sign of 
dispersion and nonlinear terms. 

In full accordance with the theoretical predictions the bright solitons have been observed 
for nonlinear backward volume waves case@ -i (in-plane field is directed parallel to the carrier 
wave vector and propagation velocity of envelope soliton), while the dark solitons are created 
in case of nolinear surface wavesll-0 (carrier wave vector and group velocity is perpendicular 
to the magnetic field). It should be especially noted that all the mentioned solitons are 
observed in narrow strips. In such geometries the transverse instabilities do not develop and 
experiments show the stable propagation of ID solitons along the waveguides. At the same 
time, in wide samples ID solitons are in general unstable0,0 and form metastable spin wave 
bullets which decay either after edge reflection or mutual interactionllljlil. 

We emphasize that the solitons in in-plane magnetized films are studied both theoreti- 
cally and experimentally only for two particular cases when the pulse propagates along or 
perpendicular to the magnetic field. Only very recently the general case of linear and non- 
linear magnetostatic wave propagation in wide samples was investigated^ for a wide range 
of angles between propagation velocity and magnetic field. In this connection the natural 
question arises: why one does not consider the nonlinear pulses in waveguides which are 
not either parallel or perpendicular to in-plane static magnetic field. As we show below for 
each magnitude of internal magnetic field it is possible to choose the direction (besides 
or 90 degrees) of waveguide respect to magnetic field direction for which stable propagation 
of envelope solitons is allowed (see inset of Fig. 1 for a geometry of the problem). We de- 
termine the limits for magnitude of magnetic field, angle between waveguide and magnetic 
field vector and pulse frequency necessary for creation of dark or bright envelope solitons. 
We also calculate their widths and propagation velocities and claim that such solitons could 
be experimentally observed. 
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II. BASIC CONSIDERATION: LINEAR MAGNETOSTATIC WAVES 



The linear consideration is based upon the Damon-Eshbach formulationEI and its general- 
ization by Hurben and Pattoni for the case of arbitrary angles between a wave vector k and 
static internal magnetic field Hq. Further we will examine only so-called " near- uniform" 
case {kd <C 1, d stands for a film thickness) and derive the dispersion expansion over the 
parameter kd up to a second order. Therefore we present here only the steps necessary for 
this purpose. Consideration of the mentioned wave number range sufficiently simplifies cal- 
culations and, besides that, most of the experiments on the magnetostatic envelope solitons 
are made having such carrier wave numbers. 

Examining an in-plane magnetized ferromagnetic film with unpinned surface spins let us 
make the following definitions: z is a direction of internal static magnetic field; f indicates 
the radius vector lying in the sample plane (y, z) and x is a coordinate along the direc- 
tion perpendicular to the film plane. Then one can write down the Landau-Lifshitz and 
magnetostatic equations: 

= —g[M x H] div(H + 4vrM) = rotH = 0. (1) 

Here g is a modulus of the gyromagnetic ratio for electrons; M is a magnetization density 
vector and H is an internal total magnetic field. Introducing dynamical dimensionless quan- 
tities rh = (M — M Q )/M Q and h = (H — H Q )/H the following equations are obtained in the 
linear limit over \m\: 

dm x dm y 

= u H {m y - h y ); = -u H (m x - h x ); 

Odd 

— (uj H h x + LO M m x ) + — (uj H h y + uj M m y ) + — {u H h z ) = 0, (2) 

where ujh = gHo and ujm = ^gMo. Defining h = gradQ and searching for the solution of 
Eq. (Q) in periodical form over t (i.e. $ and m X}y being proportional to exp(—iu>t)) we get: 
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Thus from (|3p we can write down a linear solution of ([I]) in the form@: 

$ = (Ae KX + Be- KX ) e"*^*^ for |x| < |; 
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where A, B, C D are arbitrary constants at the present stage 
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and let us remind that two dimensional vectors r— (y, z) and k = (k y , k z ) lye in the film 
plane and k = ^k* + k\. 

If k 2 > we deal with a so-called surface mode, otherwise k = i\J — k 2 and a volume mode 
exists. However, note that due to the condition of "near- uniformity" W<1 the difference 
between these two modes is negligible. 

The dispersion relation can be obtained from (|4j) if we remember about the boundary 
conditions. Particularly the functions $ and h x + 4irm x should be continuous on the bound- 
aries — d/2 and d/2 of the film. The dispersion relation for both modes could be written as 
follows: 

2k( X i + 1)k . - k 2 vX l + k 2 + k 2 ( Xi + l) 2 = 0. (5) 



Working in the limit kd <C 1 and keeping only the terms up to the second order of this 
parameter the following expression is obtained: 
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where ojq = u(k = 0) = Jujh(ojh + wm)- Then we get from @ the following expressions for 



the derivatives of oj over k y and k z : 
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(7) 



The higher approximation terms Oj(kd) and O'j(kd) are not presented here because of their 
rather cumbrous form, but we use them in the calculation as far as the leading terms in 
expressions (0) vanish in vicinity of some points, e.g. k y = or k z = 0. 
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III. WEAKLY NONLINEAR LIMIT: SOLITON SOLUTIONS 



Defining a wave envelope u 



-i(iot—kr) 



m x + im y — u • e 

and following the well known modulation approach!!!,! the nonlinear equation for wave en- 
velope u is derived (we redirect reader for details of obtaining this equation to the recent 
paper Ref.EH where the full procedure is well described): 
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(8) 



All of the coefficients are defined by formulas (|?|) except the nonlinear coefficient iV which 
could be easily calculated taking into account that in nonlinear case we have a following 
identity: 

u M = 4ngM m z . 

Then Substituting in the expansion of m z in weakly nonlinear limit m z = 1 — \u\ 2 /2 and 
using expression for N from Refs.0,i,@ we get in large wavelength limit (kd <C 1) 
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(9) 



Let us mention that if carrier wave vector is parallel or perpendicular to the static internal 
magnetic field the coefficients v y and w yz are equal to zero!,!. But for arbitrary angles 
between k and H Q that is not the case. Therefore we should introduce a new frame of 
references in order to vanish the nondiagonal term with coefficient uo'y Z . This could be done 
rotating the frame of references yz by the angle d 



£ = z cos i? + y sin d, rj = y cos i? — z sin 



where 
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where 



R = u" zz cos 2r d + 2u" zy cos i d sin -d + u! yy sin 2, d, 
S = u zz sin 2 {} — 2uj" y cos , d sin i? + u yy cos 2, d, 



v z cos$ + VySind, 
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(14) 



Afterwards in the moving frame of references 

£ 1 =£-Vit 7)1=7}- v 2 t 



5 



we come to the 2D nonlinear Schrodinger equation: 



.du I B 2 u 1 „d 2 u „ rl l9 







(15) 



and can write down its ID bright or dark soliton solutions assuming that soliton envelope 
is a function only of variables £i and t (thus soliton propagates along a spatial axis £ with 
a velocity vi). If iVi? < we have bright soliton with envelope 
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while in case NR > dark soliton solution is permitted: 
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where A denotes the contrast of dark soliton (if A = 1 one has a black dark soliton and gray 
dark otherwise) and soliton width A is defined for both cases by the same way: 
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(18) 



Now we shall discuss the question about the stability of these ID solitons. 



IV. STABLE SOLITONS IN WAVEGUIDES FOR TILTED MAGNETIC FIELDS 

As well known ID soliton solutions ( ]TB| ) and ( |I7| ) of 2D NLS are not stable to the trans- 
verse modulations with wavenumbers < k < K c . According to the recent results (see e.g. 
Ref.0) k c ~ 1/A, thus, if the limits of transverse variable rji is less than soliton width the 
instabilities do not develop and ID soliton solutions (IB) and flT7p would be stable. When 



one has fully spatial transverse variable r)i the above condition means that narrow samples 
should be used. In our case we have the mixed variable r]i = r] — v and therefore also the 
condition for time dependent part has to be introduced: v 2 t < A and in case 

v 2 = (19) 

the transverse instabilities do not develop even for infinite time. Thus besides the condition 
(|TT1) we get from (0) and (|19D additional condition on the stable soliton parameters: 



tgti = X (20) 



Afterwards, in view of both conditions (0) and (E^) we finally obtain the following equality: 



Solving (|2lD as an expansion over small parameter kd we simply come to the following 
expression for the angle <p between carrier wavevector k and static magnetic field: 

sm(p=—A H 5-— fcd . (22) 

V u H + u M \ 3(lu 2 h - uj 2 m ) J 

Obviously there exist also trivial solutions ip = and ip = 90° which will not be considered 
as long as they correspond to the well studied cases of bright backward volume wave and 
dark surface wave solitons, respectively. Further using the condition ( pH|) and definitions 
for dispersion coefficient (|1^) we get the expressions for •& and R as functions of expansion 
parameter kd: 

sh*=J^»—(l ^—kd); R= * (23) 

V U) M + uj h V UJ H + UJ M J 2^uj h {uj m + U) H ) k 

As long as the nonlinear coefficient N according to @ is always negative the possibility of 
appearance of dark or bright solitons depends on the sign of dispersion coefficient R. In view 
of the second relation in ( p3|) we can conclude that bright solitons appears if ujh < ujm while 
dark solitons could be created for larger magnetic fields wg > ujm- From Exps. ([Tj|) and 
(|) we are also able to get the expressions for soliton width and detuning of pulse frequency, 
respectively: 



d 2\uj m — uj h \ % / uJm ujm n ,x 2 fnA \ 
A= M V vT A — 5 w-w = — J {kd) , (24) 

\U\max V UJ H kd 3 \ UJ H + UJ M UJ H ~ UJ M 

while the soliton propagation velocity could be given by simple approximate formula: 



uJAid / ujm /f ,_s 

t> = ~ , / 25 

2 V^ + Wh V ; 

It should be noted that our perturbative approach violates if uj h — > ujm- Besides that, we 
have a following restriction on the internal static magnetic field: ujh > 0.3ujm- Otherwise 
the threshold for three magnon processes is reached and localized nonlinear wave will decay 
rapidlyi. 

As we see all of the quantities u — u> and A specifying the soliton are the functions 
of kd and h = uj h /ujm- Thus it is possible to plot uj — uj and A as the functions of •& for 
different h (see Figs. 1 and 2). 

In Fig. 1 we present how to choose the sample geometry (in other words how to choose 
an angle 1? between waveguide and static magnetic field) and frequency of applied pulse for 
various h = ujh/ujm in order to create bright or dark soliton. While in Fig. 2 we show the 
dependence of soliton width on the geometry of the problem and static magnetic field. In 
both cases the curves are limited because of restriction of "near uniformity" kd <S 1 and the 
following parameters for YIG film are used: d = 10/im and ujm = 1750 Oe. Note that bright 
solitons appear when the angle between waveguide and static magnetic field is less than 45°. 
Besides that, the detuning of pulse should be positive in order to create dark solitons. For 
the purpose to create bright solitons the angles should be larger than 45° and detuning has 
to be negative. 
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V. CONCLUSIONS AND POSSIBLE EXPERIMENTAL SETUP 



Summarizing we can declare that the possibility of stable magnetostatic soliton prop- 
agation in in-plane magnetized ferromagnetic films in presence of tilted (from waveguide 
direction) static in-plane magnetic fields is proved. The widths and velocities as well as 
the range of angles between waveguide and magnetic field is obtained for which the stable 
soliton propagation is allowed. 

However, As it was mentioned by the anonymous referee the waveguide border (see dashed 
line in the inset of Fig. 1) could cause the reflection of carrier wave (wave vector k) what 
will change the group velocity destructing thus the soliton. To avoid such a possibility we 
propose to use tube like magnetic waveguides (see Fig. 3). Then the carrier wave will not 
be reflected and, besides that, the condition of quasi-one dimensionality still holds. Let us 
make the following choice of the parameters of the problem: tube diameter L = 0.5mm; film 
thickness d = 10/xm and 1/L <C k -C 1/d. Thus the near uniformity condition kd <C 1 is still 
valid and simultaneously the consideration of the sample as locally flat is allowed proving 
thus approximate validity of solutions like Exp. (£|). 
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FIG. 1: Detuning of the pulse frequency oj — uq versus an anjjfe [^Oegfe^sJ^aveguide and static 
in-plane magnetic field for its different magnitude h = ljh/ujm- Dashed and solid lines correspond 
to the dark and bright soliton cases, respectively. Inset shows the geometry of the problem and 
dashed lines indicate the direction of a waveguide. The following film parameters are used in 
calculations: Film thickness d = 10/um and value of demagnetizing field Hm = 1750 Oe. 
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FIG. 2: Dependence of the envelope soliton width on the angle 6 between waveguide direction 
and in-plane magnetic field for its different value. Dashed and solid lines indicate dark and bright 
soliton cases, respectively. As in the previous figure h = uh/wm\ Hm = 1750 Oe; film thickness is 
equal to d = 10/xm and, besides that, relative amplitude of soliton is taken as follows: |n max | = 0.1. 
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FIG. 3: The possible experimental setup for observation of stable magnetostatic solitons. In each 
point of the sample magnetic field is parallel to the tube surface and tilted by the angle 6 from the 
tube symmetry axis which is parallel to the envelope soliton propagation direction. 
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